R Notebook

library(dplyr)

Registered S3 method overwritten by 'dplyr':
  method           from
  print.rowwise_df     

Attaching package: ‘dplyr’

The following objects are masked from ‘package:stats’:

    filter, lag

The following objects are masked from ‘package:base’:

    intersect, setdiff, setequal, union

library(devtools)

Loading required package: usethis

library(leaflet)

Registered S3 methods overwritten by 'htmltools':
  method               from         
  print.html           tools:rstudio
  print.shiny.tag      tools:rstudio
  print.shiny.tag.list tools:rstudio
Registered S3 method overwritten by 'htmlwidgets':
  method           from         
  print.htmlwidget tools:rstudio

1. Location of Severe Fires

building_fire <- read.csv("/Users/carolineyu/Documents/1111/course_materials/Exercises/07_fire/building_fires.csv")
FDNY <- read.csv("/Users/carolineyu/Documents/1111/course_materials/Exercises/07_fire/FDNY_Firehouse_Listing.csv")

building_fire$HIGHEST_LEVEL_DESC <- case_when( building_fire$HIGHEST_LEVEL_DESC == "0 - Initial alarm" ~"0 - Initial alarm",
building_fire$HIGHEST_LEVEL_DESC == " 11 - First Alarm " ~ "1 - Initial alarm, less than 10-75",
building_fire$HIGHEST_LEVEL_DESC == " 1 - More than initial alarm, less than Signal 7-5" ~ " 1 - More than initial alarm, less than Signal 7-5 ",
building_fire$HIGHEST_LEVEL_DESC %in% c("2 - 2nd alarm ","22 - Second Alarm") ~ "2 - 2nd alarm",
building_fire$HIGHEST_LEVEL_DESC %in% c("3 - 3rd alarm ", " 33 - Third Alarm ") ~ "3 - 3rd alarm",
building_fire$HIGHEST_LEVEL_DESC%in% c("4 - 4th alarm", " 44 - Fourth Alarm") ~ "4 - 4th alarm",
building_fire$HIGHEST_LEVEL_DESC %in% c(" 5 - 5th alarm", " 55 - Fifth Alarm ") ~ " 5 - 5th alarm",
building_fire$HIGHEST_LEVEL_DESC %in% c("7 - Signal 7-5", "75 - All Hands Working") ~ "7 - Signal 7-5")
Highest_severity <- building_fire %>%
  filter(HIGHEST_LEVEL_DESC == "7 - Signal 7-5")

content <- paste("Severity:", Highest_severity$HIGHEST_LEVEL_DESC,"<br/c>",
                 "When:",Highest_severity$INCIDENT_DATE_TIME,"<br/c>",
                 "Where:",Highest_severity$PROPERTY_USE_DESC,"<br/c")

m <- leaflet(Highest_severity) %>%
  addTiles()
m1 <- m %>%
  addCircles(lng = ~lon, lat = ~lat, popup = content)

2. Layers and Clusters

1. Color by Type of Property

Highest_severity$PROPERTY_USE_DESC <- case_when(grepl("dwelling",Highest_severity$PROPERTY_USE_DESC)== TRUE | Highest_severity$PROPERTY_USE_DESC %in% c("400 - Residential, other  ","439 - Boarding/rooming house, residential hotels "," 449 - Hotel/motel, commercial"," 460 - Dormitory-type residence, other ","464 - Barracks, dormitory ") ~ "Dwelling, Residential",Highest_severity$PROPERTY_USE_DESC%in% c(" 161 - Restaurant or cafeteria") ~"Restaurant", Highest_severity$PROPERTY_USE_DESC%in% c("881 - Parking garage, (detached residential garage)","965 - Vehicle parking area "," 882 - Parking garage, general vehicle  ") ~"Parking", Highest_severity$PROPERTY_USE_DESC%in% c("564 - Laundry, dry cleaning ","579 - Motor vehicle or boat sales, services, repair","559 - Recreational, hobby, home repair sales, pet store ","500 - Mercantile, business, other   ","549 - Specialty shop"," 519 - Food and beverage sales, grocery store ","580 - General retail, other","557 - Personal service, including barber & beauty shops ","511 - Convenience store    ","539 - Household goods, sales, repairs ","183 - Movie theater ","180 - Studio/theater, other ","160 - Eating, drinking places, other ","529 - Textile, wearing apparel sales","581 - Department or discount store","186 - Film/movie production studio","181 - Live performance theater","254 - Day care, in commercial property","161 - Restaurant or cafeteria","142 - Clubhouse ","140 - Clubs, other","143 - Yacht Club ","162 - Bar or nightclub","144 - Casino, gambling clubs","141 - Athletic/health club") ~"Commercial Property", Highest_severity$PROPERTY_USE_DESC%in% c(" 131 - Church, mosque, synagogue, temple, chapel")  ~ "Religious Places", Highest_severity$PROPERTY_USE_DESC%in% c(" 173 - Bus station     ","170 - Passenger terminal, other ","174 - Rapid transit station ")  ~ "Public Transportation", grepl("Clinic",Highest_severity$PROPERTY_USE_DESC) == TRUE | Highest_severity$PROPERTY_USE_DESC%in% c("331 - Hospital - medical or psychiatric ","311 - 24-hour care Nursing homes, 4 or more persons"," 332 - Hospices "," 300 - Health care, detention, & correction, other ","342 - Doctor, dentist or oral surgeon office"," 321 - Mental retardation/development disability facility"," 340 - Clinics, doctors offices, hemodialysis cntr, other","363 - Reformatory, juvenile detention center") ~ "Clinic",
Highest_severity$PROPERTY_USE_DESC%in% c("210 - Schools, non-adult, other","241 - Adult education center, college classroom  ","215 - High school/junior high school/middle school","211 - Preschool","213 - Elementary school, including kindergarten")  ~ "School", Highest_severity$PROPERTY_USE_DESC%in% c("123 - Stadium, arena"," 931 - Open land or field ","962 - Residential street, road or residential driveway  ","322 - Alcohol or substance abuse recovery center","610 - Energy production plant, other","926 - Outbuilding, protective shelter ","000 - Property Use, other"," 981 - Construction site","459 - Residential board and care  "," 152 - Museum ","130 - Places of worship, funeral parlors, other","963 - Street or road in commercial area ","648 - Sanitation utility","110 - Fixed-use recreation places, other","571 - Service station, gas station ","182 - Auditorium, concert hall  "," 808 - Outbuilding or shed "," 121 - Ballroom, gymnasium  ","112 - Billiard center, pool hall ","569 - Professional supplies, services ","124 - Playground","635 - Computer center","134 - Funeral parlor","984 - Industrial plant yard - area","629 - Laboratory or science lababoratory","936 - Vacant lot","250 - Day care, other (Conversion only)","615 - Electric-generating plant ","639 - Communications center","898 - Dock, marina, pier, wharf ","952 - Railroad yard ")  ~ "Public Area",Highest_severity$PROPERTY_USE_DESC%in% c("891 - Warehouse","899 - Residential or self-storage units ","880 - Vehicle storage, other  ","800 - Storage, other","839 - Refrigerated storage ","807 - Outside material storage area")  ~ "Warehouse/Storage",Highest_severity$PROPERTY_USE_DESC%in% c(" 599 - Business office","593 - Office:  veterinary or research ","150 - Public or government, other   ","596 - Post office or mailing firms","365 - Police station","155 - Courthouse","Bank")  ~ "Office")

pal = colorFactor("Set1", domain = Highest_severity$PROPERTY_USE_DESC)
color_pro = pal(Highest_severity$PROPERTY_USE_DESC)
m2 <- m %>%
  addCircles(lng = ~lon, lat = ~lat, color = color_pro,popup = ~as.character(Highest_severity$PROPERTY_USE_DESC)) %>%
addLegend(pal = pal, values = ~ Highest_severity$PROPERTY_USE_DESC, title = "Type of Property")
m2

2. Cluster

m3 <- m %>%
addCircleMarkers(lng = ~lon, lat = ~lat,color = color_pro, popup = content, clusterOptions = markerClusterOptions())
m3

3. Fire House

m_l1 <- leaflet() %>%
addTiles()%>%
addCircleMarkers(data = Highest_severity ,group = "Incidents", radius  = ~ Highest_severity$UNITS_ONSCENE,popup = content) %>%
addCircleMarkers(data = FDNY,group = "Firehouses",color = "#F05",popup = ~as.character(FDNY$FacilityName),clusterOptions = markerClusterOptions()) %>%
addLayersControl(
baseGroups = c("NYCmap"),
overlayGroups = c("Incidents","Firehouses"))

Assuming "lon" and "lat" are longitude and latitude, respectively
Assuming "Longitude" and "Latitude" are longitude and latitude, respectively
Data contains 5 rows with either missing or invalid lat/lon values and will be ignored

m_l1

4. Distance from Firehouse and Response Time

1. Calculate Distance

building_fire1<- building_fire %>%
mutate(distime = difftime(strptime(as.character(building_fire$ARRIVAL_DATE_TIME),"%m/%d/%Y %I:%M:%S %p"),strptime(as.character(building_fire$INCIDENT_DATE_TIME),"%m/%d/%Y %I:%M:%S %p"),units = "secs"))

library(geosphere)
library(SoDA)
library(ggplot2)

fire_lon_lat <- data.frame(building_fire1$lon, building_fire1$lat)
FDNY_lon_lat <- data.frame(FDNY$Longitude, FDNY$Latitude)
a <-  distm(as.vector(fire_lon_lat),as.vector(FDNY_lon_lat),fun = distGeo)

building_fire2 <- building_fire1 %>%
mutate(mindist = NA)
for (i in 1:nrow(a)) {
building_fire2$mindist[i] <- min(a[i,],na.rm = TRUE)
}

p_distance <- ggplot(building_fire2, aes(x = mindist, y = distime)) + geom_point() + xlim(0,5000) + ylim(0,1000) + facet_wrap(as.factor(building_fire2$HIGHEST_LEVEL_DESC)) + geom_smooth() + labs(x="Nearest Firehouse in meters", y = "First arrival in secs", title = "Distance From Firehouse V.S. Response Time")
p_distance

2. Map of Response Times

library(rgdal)

Loading required package: sp
rgdal: version: 1.4-8, (SVN revision 845)
 Geospatial Data Abstraction Library extensions to R successfully loaded
 Loaded GDAL runtime: GDAL 2.4.2, released 2019/06/28
 Path to GDAL shared files: /Library/Frameworks/R.framework/Versions/3.6/Resources/library/rgdal/gdal
 GDAL binary built with GEOS: FALSE 
 Loaded PROJ.4 runtime: Rel. 5.2.0, September 15th, 2018, [PJ_VERSION: 520]
 Path to PROJ.4 shared files: /Library/Frameworks/R.framework/Versions/3.6/Resources/library/rgdal/proj
 Linking to sp version: 1.3-2 

library(geojson)

Attaching package: ‘geojson’

The following object is masked from ‘package:graphics’:

    polygon

library(geojsonio)

Attaching package: ‘geojsonio’

The following object is masked from ‘package:geosphere’:

    centroid

The following object is masked from ‘package:devtools’:

    lint

The following object is masked from ‘package:base’:

    pretty

nycity <- readOGR("/Users/carolineyu/Desktop/borough_boundaries.geojson", verbose = FALSE)

content_time <- paste("Severity:", building_fire2$HIGHEST_LEVEL_DESC,"<br/c>",
"When:",building_fire2$INCIDENT_DATE_TIME,"<br/c>",
"Property Type:",building_fire2$PROPERTY_USE_DESC,"<br/c",
"Response Time:",building_fire2$distime,"<br/c")

building_fire3 <- building_fire2 %>%
  group_by(BOROUGH_DESC, HIGHEST_LEVEL_DESC)%>%
  summarise(avg = mean(distime, na.rm = TRUE))

building_fire3$BOROUGH_DESC <- substring(building_fire3$BOROUGH_DESC, 5)

nycity_severity <- left_join(nycity@data,building_fire3, by = c("boro_name"="BOROUGH_DESC"))

Column `boro_name`/`BOROUGH_DESC` joining factor and character vector, coercing into character vector

mypalette <- colorNumeric(palette ="YlOrRd", domain = as.numeric(building_fire3$avg),c(180,250))

NYC_Fires <- leaflet() %>% 
  addTiles() %>%
  addPolygons(data = nycity,fillOpacity = 1,fillColor = ~mypalette(as.numeric(building_fire3$avg)))
NYC_Fires

library(lubridate)
building_fire2$YEAR <- substring(building_fire2$INCIDENT_DATE_TIME,7,10)

yearly <- building_fire2 %>%
  group_by(YEAR, BOROUGH_DESC) %>%
  summarise(avg = mean(distime, na.rm = TRUE))

library(reshape )
yearly$BOROUGH_DESC <- substring(yearly$BOROUGH_DESC, 5)

nyc_by_year <- full_join(yearly, nycity@data, by = c("BOROUGH_DESC" = "boro_name" ))

mypalette1 <- colorNumeric("YlOrRd", domain = as.numeric(nyc_by_year$avg), c(180,250))

year2013 <- left_join(nycity@data, yearly[yearly$YEAR == 2013, ], by = c("boro_name" = "BOROUGH_DESC"))

Column `boro_name`/`BOROUGH_DESC` joining factor and character vector, coercing into character vector

nyc2013 <- leaflet()%>%
  addTiles() %>%
  addPolygons(data = nycity, fillColor = ~ mypalette1(as.numeric(year2013$avg)),
              fillOpacity = 1)%>%
  addLegend(pal = mypalette1, values = as.numeric(nyc_by_year$avg), title = "YEAR 2013", "topright")

year2014 <- left_join(nycity@data, yearly[yearly$YEAR == 2014, ], by = c("boro_name" = "BOROUGH_DESC"))

Column `boro_name`/`BOROUGH_DESC` joining factor and character vector, coercing into character vector

nyc2014 <- leaflet()%>%
  addTiles() %>%
  addPolygons(data = nycity, fillColor = ~ mypalette1(as.numeric(year2014$avg)),
              fillOpacity = 1)%>%
  addLegend(pal = mypalette1, values = as.numeric(nyc_by_year$avg), title = "YEAR 2014", "topright")

year2015 <- left_join(nycity@data, yearly[yearly$YEAR == 2015, ], by = c("boro_name" = "BOROUGH_DESC"))

Column `boro_name`/`BOROUGH_DESC` joining factor and character vector, coercing into character vector

nyc2015<- leaflet()%>%
  addTiles() %>%
  addPolygons(data = nycity, fillColor = ~ mypalette1(as.numeric(year2015$avg)),
              fillOpacity = 1)%>%
  addLegend(pal = mypalette1, values = as.numeric(nyc_by_year$avg), title = "YEAR 2015", "topright")

year2016 <- left_join(nycity@data, yearly[yearly$YEAR == 2016, ], by = c("boro_name" = "BOROUGH_DESC"))

Column `boro_name`/`BOROUGH_DESC` joining factor and character vector, coercing into character vector

nyc2016 <- leaflet()%>%
  addTiles() %>%
  addPolygons(data = nycity, fillColor = ~ mypalette1(as.numeric(year2016$avg)),
              fillOpacity = 1)%>%
  addLegend(pal = mypalette1, values = as.numeric(nyc_by_year$avg), title = "YEAR 2016", "topright")

year2017 <- left_join(nycity@data, yearly[yearly$YEAR == 2017, ], by = c("boro_name" = "BOROUGH_DESC"))

Column `boro_name`/`BOROUGH_DESC` joining factor and character vector, coercing into character vector

nyc2017 <- leaflet()%>%
  addTiles() %>%
  addPolygons(data = nycity, fillColor = ~ mypalette1(as.numeric(year2017$avg)),
              fillOpacity = 1)%>%
  addLegend(pal = mypalette1, values = as.numeric(nyc_by_year$avg), title = "YEAR 2017", "topright")

year2018 <- left_join(nycity@data, yearly[yearly$YEAR == 2018, ], by = c("boro_name" = "BOROUGH_DESC"))

Column `boro_name`/`BOROUGH_DESC` joining factor and character vector, coercing into character vector

nyc2018 <- leaflet()%>%
  addTiles() %>%
  addPolygons(data = nycity, fillColor = ~ mypalette1(as.numeric(year2018$avg)),
              fillOpacity = 1)%>%
  addLegend(pal = mypalette1, values = as.numeric(nyc_by_year$avg), title = "YEAR 2018", "topright")

library(mapview)
leafsync::sync(nyc2013,nyc2014,nyc2015,nyc2016,nyc2017,nyc2018)